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Abstract 

The NA60 experiment at the CERN SPS has studied <j> meson production in In-In collisions at 158A GeV 

via both the K + K~ and the fj, + [i~ decay channels. The yields and inverse slope parameters of the mj 

spectra observed in the two channels are compatible within errors, different from the large discrepancies 

seen in Pb-Pb collisions between the hadronic (NA49) and dimuon (NA50) decay channels. Possible physics 
implications are discussed. 

Keywords: Heavy Ion Collisions, (f> puzzle 



The theory of Quantum Chromo-Dynamics (QCD) 
predicts that matter under extreme conditions of 
temperature and energy density undergoes a phase 
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transition from hadronic matter to a plasma of de- 
confined quarks and gluons (QGP). The occurrence of 
this phase can be studied in the laboratory by means 
of high energy heavy-ion collisions. Strangeness en- 
hancement has been suggested as a signature of QGP 
formation [lj]. The 4> meson, being composed of an 
ss pair, is an ideal probe for the study of strangeness 
production. Moreover, the </> mass and branching ra- 
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Table 1: Ccntrality bins, measured number of <j>, ratio between background and signal, measured <j> mass and width. 



dN c h/drj range 


{Npart) 


A^,(meas.) 


B/S 


m (MeV/c*) 


a m (MeV/c^) 


29-42 
42-95 


41 
78 


1.5 • 10 4 
1.3 ■ 10 5 


188 
240 


1019.9 ± 1.7 


6.5 ± 1.2 


96-160 


133 


6.5 • 10 5 


310 


1019.9 ± 1.2 


7.5 ± 1.1 


161-250 


177 


6.8 ■ 10 5 


465 


1019.3 ±1.8 


7.4 ±2.3 



Table 2: (p inverse slope and yield extracted from the analyses in the muon and kaon channels as a function of centrality. 



{N part ) (MeV) (0) A1A1 ( PT >O.9GeV) T C K S K (McV) (0}kk{pt > 0-9 GcV) {&kk{pt > 0) 



15 


209 ±4 


0.044 ±0.002 ±0.005 








41 


232 ±4 


0.197 ±0.009 ±0.021 




0.16 ±0.05 ±0.08 


0.64 ±0.17 ±0.03 


78 


245 ±4 


0.48 ±0.02 ±0.05 




0.47 ±0.06 ±0.02 


1.60 ±0.18 ±0.06 


133 


250 ±4 


1.03 ±0.03 ±0.07 


253 ±11 ±5 


1.01 ±0.07 ±0.05 


3.5 ±0.2 ±0.2 


177 


249 ±5 


1.65 ±0.06 ±0.07 


254 ± 13 ±6 


1.49 ±0.11 ±0.07 


4.6 ±0.3 ±0.4 



tios for the decay into kaon and lepton pairs may be 
modified in the medium 0, 0] ■ 

4> meson production was first studied at the CERN 
SPS by the NA49 [li] and NA50 0,0 experiments. 
NA49 detected the </> through its decay into K + K~ 
pairs, while NA50 studied the 4> —> t 1+ ^~ channel. 
The results obtained by the two experiments in Pb- 
Pb collisions show discrepancies both in the yield and 
in the inverse slope parameter T e g of the wit spectra. 
The (f> multiplicity in central Pb-Pb collisions mea- 
sured by NA50 in the dimuon channel is higher by 
about a factor of 4 with respect to the corresponding 
NA49 measurement in the K + K~ channel. The T c g 
value found by NA50 is about 220-230 MeV, showing 
a mild dependence on centrality, while NA49 mea- 
sured an inverse mx slope that increases with cen- 
trality and saturates at T e g ~ 300 MeV. It has to be 
stressed that the NA50 acceptance is limited to high 
Pt, while NA49 is dominated by low pt- 
The discrepancy between the NA49 and NA50 re- 
sults, known as "</> puzzle", triggered a considerable 
effort to explain the observed differences [1,0, [13, U|- 
It was argued that in-medium effects may affect the 
spectral function of the <p, causing a modification of 
its mass and partial decay widths. Moreover, kaon 
absorption and rescattering in the medium can re- 
sult in a loss of signal in the region of the 4> invariant 
mass in the K + K~ channel, thus reducing the ob- 
served yield. This effect would be concentrated at 
low pt, causing a hardening of the px spectrum in 



this channel 11, 9]. Nevertheless, according to those 
calculations, the yield in lepton pairs is expected to 
exceed the one in kaon pairs by about 50%, which 
is much lower than the observed differences. More 
recently, the CERES experiment at the SPS studied 
4> production in Pb-Pb collisions both in the K + K~ 
and dielectron channels, finding an agreement with 
the NA49 results [Hj]. However, CERES' measure- 
ment of the T e g parameter is affected by a large sta- 
tistical error that does not allow to draw firm conclu- 
sions. 

The NA60 experiment measured <f> production in In- 
In collisions at the CERN SPS. The <fi meson is de- 
tected through its decay both in muon and kaon pairs. 
In this paper, results on the 4> — > K + K~ channel will 
be presented and compared to the already published 
ones for the </> — > n + [iT decay mode [IH and to the 
existing SPS measurements in different systems. 
The detectors relevant for the present work are a 
muon spectrometer (MS) inherited from NA50 and 
a vertex telescope (VT). They are fully described 
in 3, IH • The muon spectrometer is composed of a 
toroidal magnet, 8 multi-wire proportional chambers 
for the tracking and a set of scintillator hodoscopes 
that provides the main trigger, which selects muon 
pairs coming from a common vertex. It is preceded 
by a 12 A/ thick hadron absorber. The VT is a 
high-granularity, radiation tolerant Si pixel detector, 
placed between the target and the absorber in a 2.5 T 
dipole field. It is used to determine the primary ver- 
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tex, the charged particle multiplicity and to measure 
the momentum of the charged tracks. 
The sample used for the results presented in this pa- 
per was collected with a 158 A-GeV In beam im- 
pinging on a 0.17 A/ thick In target, composed of 7 
subtargets placed in vacuum. The acquired statis- 
tics consists of 230 million triggers, mainly dimuons. 
More than 99% of the dimuon events consist of un- 
correlated pairs coming from the decay of pions and 
kaons in muons, or non-muon tracks. Since the tracks 
used in the present analysis are reconstructed using 
the Si tracker alone, the dimuon trigger acts like a 
minimum bias trigger. The distortion of the multi- 
plicity distribution introduced by the dimuon trigger- 
drops out, since the results are obtained normalizing 
the number of <fi — > K + K~ obtained in each multi- 
plicity class to the number of events in that class, as 
discussed later in the text. Therefore, no bias is in- 
troduced by the dimuon trigger. 
The charged tracks associated to the vertices recon- 
structed in the vertex tracker are used for this analy- 
sis. In order to avoid events with reinteractions of nu- 
clear fragments in the subsequent targets, only events 
with one vertex in the target region are selected. The 
vertices are required to lay inside one of the Indium 
subtargets. In addition, this cut eliminates events 
with ambiguous vertices and pile-up events, and re- 
jects about 55% of the statistics. 
The tracks are selected requiring that the reduced 
X 2 of the track fit is lower than 3. The collision 
ccntrality is determined through the measurement 
of the charged particle multiplicity. The number 
of produced charged particles N c h is extracted from 
the raw charged track multiplicity, applying a cor- 
rection that takes into account the detector accep- 
tance, reconstruction efficiency and secondary par- 
ticle production. The number of participants is 
then obtained assuming N vart oc dN c h/drj. We find 



N 



part 



dN c h/dr]\3.r |13j], with a systematic error of 



about 10% for peripheral collisions and 5% for cen- 
tral collisions. 

The analysis in the K + K~ channel is performed as- 
suming that all the charged particles associated with 
the primary vertex are kaons, and building all the 
possible opposite sign pairs among the tracks of each 
event (as discussed below, like sign pairs are used 



for systematic checks and normalization purposes). 
This results in a considerable combinatorial back- 
ground, which is reduced by a factor 10 in the <f> 
mass range applying a cut on the pair's opening an- 
gle, 0.005 < 9kk < 0.15 rad. Due to the low ac- 
ceptance for the signal at low pt, and to the cor- 
responding high background to signal ratio (about 
three times higher than at high px) a cut on the pair 
transverse momentum pt > 0.9 GeV/c was applied. 
Furthermore, in order to exclude the borders of the 
detector acceptance, only tracks in the rapidity re- 
gion 2.9 < y < 3.7 are selected, where the rapidity is 
calculated assigning the kaon mass to the tracks. Af- 
ter applying these cuts, the ratio between background 
and signal ranges from ~ 190 to ~ 460, depending on 
centrality. 

The residual background is subtracted with an event 
mixing technique: events from runs taken in ho- 
mogeneous conditions are grouped in pools accord- 
ing to the position of the target associated to the 
vertex, the centrality of the collision and the direc- 
tion of the event plane (16|. The event plane az- 



imuthal angle is calculated by means of the Bow vec- 
tor. Its components for the second harmonic are 
Q% = Si' 10 * cos(2tfi) and Q\ = X\i0jsin(2</jj) 
where the weight Wi associated to the i — th track 
is its transverse momentum and ipi is its azimuthal 
angle. Tracks from different events belonging to the 
same pool are mixed. Each event is mixed with two 
other events. Mixed events are subject to the same 
cuts as the ones applied to the real data. 
The obtained mixed spectra are normalized asking 
that in the mass region 1.02 < m < 1.06 GeV/c 2 the 
number of mixed like-sign pairs coincides with the 
corresponding one in the real data. Changing the 
mass range yields negligible differences in the results. 
Alternatively, the normalization factor is chosen such 
that the integral of the opposite sign mass spectra in 
the whole mass range is equal for the real and mixed 
samples. For both methods, the normalization fac- 
tor is evaluated either for each subtarget separately, 
or averaged over the targets. The differences in the 
results coming from the choice of the normalization 
criterion, causing variations of about 4% in the de- 
termination of the number of (f> mesons per centrality 
bin, are taken into account in the systematic error. 
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Figure 1: Invariant mass spectrum of the opposite sign pairs af- 
ter combinatorial background subtraction for pt > 0.9 GeV/c 
integrated in centrality. 



The invariant mass spectrum of the opposite 
sign pairs after combinatorial background subtrac- 
tion (normalized with the like-sign pairs) is shown in 
Fig. [I] integrated in centrality. The function for the 
<fi peak is determined fitting the mass distribution ob- 
tained with an overlay Monte Carlo simulation, that 
consists in reconstructing a generated <j) decaying into 
a kaon pair on top of a real event. We use a gaussian 
superimposed to an empirical function that takes into 
account the mass tails and accounts for about 5% 
of the total number of <fr. Several functions and fit 
mass ranges have been tested to describe the resid- 
ual background, including polynomials of first and 
second order and functions that are null at the KK 
mass threshold. In order to check that the fit func- 
tions correctly reproduce the signal and the residual 
background, the fit procedure has been applied to 
the like-sign invariant mass spectra in several Pt in- 
tervals, after the subtraction of the event-mixed spec- 
tra. The like-sign pairs spectra have a shape which is 
very similar to the corresponding opposite sign ones 
for masses above the <j> peak (m > 1.04 GeV/c 2 ), 
and show a smooth trend down to to = 1 GeV/c 2 . 
Since these pairs do not contain any signal, then a fit 



applied to them should result in a null signal com- 
ponent, while the presence of a fictitious peak would 
indicate a bias in the fitting procedure. The fit ap- 
plied to the like-sign pairs gives a signal component 
compatible with zero in all the pt bins, showing that 
no artifact is introduced by the choice of the fit func- 
tion. 

It has to be noted that if the <f> peak position and 
width are left as free parameters of the fit, the corre- 
sponding values are = 1019.5 ± 0.3 ± 1.2 MeV/c 2 
and a m = 7.8 ± 0.3 ± 1.2 MeV/c 2 , in good agree- 
ment with the simulations. Corresponding results as 
a function of centrality are reported in Table [I] to- 
gether with the number of reconstructed <j> mesons, 
obtained by integrating the function describing the 
4> peak, and the ratio between background and sig- 
nal in the mass range 1.005 < to < 1.035 GeV/c 2 . 
Statistical and systematic errors for and a m are 
added in quadrature. The values show no dependence 
on centrality, indicating that no modification of the 
spectral function is visible within the sensitivity of 
the measurement. A similar result is obtained in the 
dimuon channel 13]. The results shown in the follow- 
ing are obtained fixing to^ and a m to 1.019 GeV/c 2 
and 7.8 MeV/c 2 , according to the Monte Carlo val- 
ues. 

In order to extract the m<r distributions, the fit to 
the invariant mass spectra is performed in several pt 
intervals having a size of 200 MeV/c. It was checked 
that results do not depend on the choice of the bin 
size. The mr distributions are then corrected for the 
geometrical acceptance and reconstruction efficiency 
with an overlay Monte Carlo simulation. The ra- 
pidity and decay angle distributions used for the 
generation are tuned to the ones measured in muons. 
The former is a gaussian with a y = 1.13, while the 
dN/d cos distribution is assumed to be flat [l3[ . 

The acceptance in multiplied by reconstruction ef- 
ficiency is almost flat for pr > 0.9 GeV/c, about 15% 
in full rapidity. A reliable tot distribution could be 
extracted only for semicentral {(N part } = 133) and 
central ((N part ) — 177) collisions, while for semipe- 
ripheral collisions ((N part ) — 78), due to statistics 
limitations, the resulting inverse slope parameter was 
not stable when varying the analysis criteria. Results 
are reported in Fig. [21 where the distributions in the 
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Figure 2: Normalized transverse mass spectra for semi-central 
(triangles) and central (circles) collisions for the cf> — > K + K~ 
channel (full symbols) compared to the corresponding ones in 
muon pairs (open symbols) |13|| . Results for central collisions 
are scaled by a factor 2 in order to improve readibility. Colour 
online. 



Figure 3: (<j>)/N par t as a function of the number of partici- 
pants in In-In collisions in the (f> — > K + K~ (full circles) and 
4> — > (open crosses) channels for > 0.9 GeV/c. Boxes: 
systematic errors. Colour online. 



4> — > fi/-j. channel in the same centrality intervals are 
reported for comparison. All the spectra are nor- 
malized to the <p multiplicity as shown in Table [5] 
(details on the method used to measure the <j> yield 
are reported below). The mj spectra are fitted with 
the function l/mTdN/drriT oc e ~ mT ^ Tatt . The re- 
duced x 2 is ~ 1 m both cases. Results are reported 
in Table [5J The systematic error is dominated by the 
choice of the function used for the background. The 
systematic errors due to the choice of the normaliza- 
tion criterion, to the variation of the analysis cuts and 
of the starting point of the fit to the mass spectra are 
also taken into account. Results are in good agree- 
ment with the ones obtained in dimuons in the full 
Pt range, also reported in Table [2j Since in presence 
of radial flow the T Q g value may depend on the px 
range, the fit to the dimuon spectra was restricted to 
the range p T > 0.9 GeV/c, giving T eff = 252 ± 4 and 
247 ± 3 MeV for semicentral and central collisions, 
still in agreement with the measurement in kaons. 

The raw </) multiplicity is determined fitting the 
mass spectra for pr > 0.9 GeV/c and dividing the 



number of <f> mesons obtained by the total number 
of events selected for this analysis. This value is 
then corrected for the branching ratio in kaon pairs, 
(49.2 ± 0.6)% [18J and for the acceptance, evaluated 
through a Monte Carlo simulation. The inverse tut 
slopes used for the calculation of the acceptance are 
the ones measured in the kaon channel, where avail- 
able; otherwise the values obtained in dimuons are 
used. Table[5]reports the results for the <fi multiplicity 
in the dimuon and kaon channels for px > 0.9 GeV/c. 
The main contributions to the systematic error are 
given by the uncertainty in the choice of the fit func- 
tion for the residual background component in the fit 
of the invariant mass spectra and the normalization 
criterion. The systematic errors due to the variation 
of the analysis cuts and of the starting point of the 
fit to the mass spectra are also taken into account. 
In Fig. [3] the ratio (<p) /N part as a function of N part 
in the <f> ~^ K + K~ channel for pt > 0.9 GeV/c is 
compared to the corresponding one in dimuons. The 
additional contribution due to the systematic error 
in N part , affecting both channels in the same way, 
is not displayed. It can be seen that the yield in 
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the hadronic channel and the one in the dileptonic 
channel are in agreement within the errors. A ratio 
between the <f> yields in dimuons and in kaons larger 
than 1.18 in the common px range is excluded at 95% 
C.L. 

The <j> multiplicity in kaons was also calculated in 
the full px range. Results are reported in the last 
column of Table [2J The systematic error due to the 
uncertainty in the extrapolation to full px caused 
by the error in the inverse slope was taken into ac- 
count. It has to be stressed that the extrapolation 
to px — was done under the hypothesis that the 
inverse slope does not change. Models like AMPT ^ 
predict that kaon rescattering causes a depletion in 
the number of reconstructed cj) — > KK which is con- 
centrated at low pt, causing an increase of T e g in 
that region. In Pb-Pb central collisions, the effect 
would be already visible for tot — to > 0.34 GeV, 
corresponding to the region covered in this analysis. 
In that region, the difference between T e ff in kaon 
and muon pairs would range from 30 to 50 MeV, and 
the fractional loss in kaons at m-r — too ~ 0.34 GeV 
would range from 35% to 50%. Such effects are not 
seen in our data, suggesting that if any rescattering 
effect is present, it would be concentrated at lower 
Pt- In the range < tot — to < 0.34 GeV the 
AMPT model predicts an inverse tot slope of about 
330 MeV, close to the value observed by NA49 in 
central Pb-Pb collisions. Since no theoretical pre- 
dictions for In-In collisions are available, we esti- 
mated the effect of the change of slope in the ex- 
trapolation to px = under the extreme hypothesis 
T cff = 330 MeV for < m T - m < 0.34 GeV, while 
for tot — too > 0.34 GeV the measured value is kept. 
This variation would cause a reduction of the <fi yield 
in kaons of about 12%. This quantity can be consid- 
ered as a conservative estimation of the uncertainty 
in the extrapolation due to a possible suppression 
mechanism of the hadronic channel. 

In order to compare to other collision systems, the 
inverse slope and the enhancement, quantified as the 
ratio (4>)/N par t in the full pt range, are plotted in 
Fig. [?] as a function of the number of participants for 
central C-C, Si-Si, In-In and Pb-Pb collisions 0,0. 
For this comparison, the systematic error on N part , of 
about 5% in central collisions, is taken into account 
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Figure 4: Inverse slope (top) and (<f>) /Npart in full px range 
(bottom) as a function of N par t for central collisions. Boxes: 
systematic errors. Colour online. 



in the calculation of the ratio ((f)) /N part for the NA60 
points. 

The inverse slope shows an initial fast increase at 
low Np ar t values, that becomes less pronounced go- 
ing towards higher N part . A lower value is observed 
by NA50 as compared both to the CERES and NA49 
measurements in Pb-Pb and to the NA60 In-In points 
in the hadronic and dileptonic channels. 

As stated above, in the presence of radial flow T e g 
depends on pt- The NA60 analysis in dimuons 13| . 
performed at low and high pt ranges, shows a dif- 
ference limited to about 15 MeV in In-In. This fact, 



complemented by the detailed blast wave analysis of 
the dimuon data discussed in makes it difficult 
to ascribe only to radial flow the large differences in 
T e ff shown in Fig. |4j A further flattening caused by 
kaon rescattering and absorption may lead to larger 
T c ff values in the hadronic channel in Pb-Pb. 
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Concerning the enhancement, the NA60 measure- 
ments in the dilepton and hadron channels can differ 
up to 18% (at 95% CL) for p T > 0.9 GeV/c and 
22% in full pt considering the uncertainty in the 
extrapolation arising from a possible suppression 
of the hadronic channel at low pr- In Fig. @] the 
NA50 result is extrapolated to full phase space 
using T eff = 220 McV, according to the NA50 
measurement in peripheral collisions. Even assuming 
as an extreme case T e g — 300 MeV, as obtained by 
NA49 in central collisions, the NA50 enhancement 
would exceed by a factor of ~ 2 the central Pb-Pb 
values measured in kaons. 

The yield per participant measured in In-In exceeds 
the one observed in Pb-Pb in the hadronic channel 
(both NA49 and CERES) by about 30%. This 
might suggest a suppressing mechanism for the kaon 
channel, below experimental sensitivity in In-In, 
which shows up in Pb-Pb. The CERES measurement 
in diclectrons is in agreement with the one in kaons 
but, given the measurement errors, it cannot rule 
out differences of the order of 40-50%, expected by 
models including kaon rescattering. If one assumes 
that no suppression mechanism affects the Pb-Pb 
collisions, it is then difficult to understand why the 
4> multiplicity in central Pb-Pb is smaller than in 
central In-In. 

In conclusion, the inverse slopes and yields measured 
in the kaon and muon channels in In-In collisions 
are in agreement, excluding for central collisions a 
difference in the yields larger than 18 % (at 95 % 
C.L.) in the common px range. In addition, no 
modification of the 4> mass and width is observed 
as a function of centrality. When the comparison is 
extended to other systems, it is difficult to reconcile 
all of the observations into a coherent picture, albeit 
there is some hint for a possible physics mechanism 
leading to a difference in the two channels in Pb-Pb 
collisions, while producing no remarkable difference 
in In-In collisions. 
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